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Abstract 

We study the large population limit of the Moran process, assuming weak- select ion, and for 
different scalings. Depending on the particular choice of scalings, we obtain a continuous 
model that may highlight the genetic-drift (neutral evolution) or natural selection; for one 
precise scaling, both effects are present. For the scalings that take the genetic-drift into 
account, the continuous model is given by a singular diffusion equation, together with two 
conservation laws that are already present at the discrete level. For scalings that take into 
account only natural selection, we obtain a hyperbolic singular equation that embeds the 
Replicator Dynamics and satisfies only one conservation law. The derivation is made in 
two steps: a formal one, where the candidate limit model is obtained, and a rigorous one, 
where convergence of the probability density is proved. Additional results on the fixation 
probabilities are also presented. 

Key words: Moran process, replicator dynamics, Kimura equations, drift-diffusion 
equations. 



1. Introduction 

1.1. Background 

An important effort in contemporary Evolutionary Biology is to understand the plethora 
of models used to describe the time evolution of a given population, here to be understood 
on a broad sense. Different models were built to model different aspects of evolution. 
Population genetics models, aiming at modeling genetic variability, had a natural start 
with discrete stochastic models at the individual level as, for instance, the Moran and 



Wright-Fisher models [12 ] . On the other hand, ecological and epidemiological models 
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assume focus on macroscopic features of the population as is the case of the SIR model 



and the Lotka-Volterra system [251 ] . This is also the case of the replicator equation [15 
and adaptive dynamics 37, 0]. Mean field approximation were also used, from the onset 



19, 12 
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Alternative modelings efforts 
graph theory 



among 



in population genetics, as in the Kimura equation 
include evolutionary models through cellular automata 
others. More recently, stochastic effects were also included at the continuous level, for 
instance in the replicator equation [l], [l7j]. Evolution of structured populations has led 
to the study of integro-differential equations 0, [3l[. A review on continuous mutation- 
selection models can be found in Models direct at the study of fixation probabilities 



can be found in 20, 21, 22, 32 



The importance of unifying evolutionary models has been already observed in 3^ . 
The difference between evolutionary game dynamics, which can lead to non-homogeneous 



distribution [15|, [34J in the long time, and the diffusive models that usually led to pure 



fixation has been also pointed out by jl, 0]. A detailed study between the relationship 
between discrete and continuous models was also carried out in 



1.2. Modeling and techniques 

The focus of this work is to understand the large population limit of discrete processes 
rigorously and to show the intimate connection between apparently diverse models. 



Thus, we consider a population 
having a fixed number iV of 



1.2.1. The discrete process 

We shall use the Moran process as a paradigm 24 
of two types — usually called the mutant and the resident 
individuals. At a given time, an individual is chosen to reproduce with probability propor- 
tional to a certain function termed the fitness, which we shall denote by The function 
\I/ is usually taken as function of the frequency of a particular type in the population, 
since this accounts for most of the biological relevant models, and it is convenient from a 
mathematical standpoint. The individual then replaces a second individual (possibly the 
same one), which is chosen with a uniform probability 1/N. Time is then increased by At, 
the time-step, and the process is repeated. 

From a theoretical point of view, while the discrete model is well understood through 
standard Markov Chain theory, the invariants of these models, other than probability 
conservation, is usually neglected, and we briefly review its mathematical properties. 

The fitness will depend, parametrically, on iV and At. When N is large, we need to 
have small time-steps in order to keep an appropriate balance. In this natural 
assumption is the so-called weak selection principle: 



*®{x,N,At) « l + (At)^«(x)+o((At) 1 



N -> oo, At -> , v > 0. 



(1) 



1.2.2. Scaling up 

We are interested in studying the so-called thermodynamical limits of the discrete 
models. While these limits can be seen as mean-field type approximation, they are more 
general than the typical diffusion approximation by partial differential equations through 
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the central limit theorem and Feymann-Kac formula. They are also somewhat less general 
than the approach by [8], for some aspects of finite population dynamics. 

More precisely, consider a one-parameter family of models Ai £ , e > 0, a certain model 
Aio and 

1. Consider a certain family of initial conditions /* such that lim e ^ fl — fl m a certain 
sense; 

2. Evolve through the model M. £ the initial condition /* and through the model M.$ the 
initial condition fl until the time t 6 (0, oo], obtaining f £ {t) and f (t) respectively, 
for t < t ; 

If, for every time t < to, we have that lim £ ^o f £ (t) = fo{t), in the same sense as before, 
then we say that the model M. £ converge in the limit e — > until time t and in the same 
sense as before to the model M.q. 

In the limit of large population, and small time-steps, we define the probability density 
of the states as: 

p(xA= lim P(zMt) (2) 

N~> oo.At^O 1 ' A 



l/N 



The function p defined by equation ([2]) will, in most cases, satisfy a certain differential 
equation, i.e., the continuous model. The particular model will be different depending on 
the precise scalings imposed. More precisely, by scaling, we mean the relationship imposed 
between N and At when passing to the limits N — * oo, At — ► 0. In this work, we also 
call "the thermodynamical limit" the limits iV — ► oo, At —>■ jointly to a given scaling. 
By extension, we also call "the thermodynamical limit" the limit model. Without loss 
of generality, we may impose scalings given by (At) 1 oc JV, n > 0. Then, we find the 
following table of possible thermodynamical limits: 



n 


V 


Large population limit 


< \x < | 




Undefined 


11 = \ 




Undefined 




" = i 


Replicator-Diffusion Equation 


M= 1 


v>\ 


Diffusion Equation 


M> 1 


< v < 1 — 


Undefined 


n>\ 


< z/ = 1 - 


Replicator Equation 


n>\ 


^ > 1 — yU 


Trivial 



This table is synthesized pictorially in Figure [U 

1.2.3. Then analysis of continuous models 

The derivation will first be presented in a formal way and then proved rigorous, once 
a candidate limit is obtained. The techniques used here come from classical numerical 
analysis [3^] rather than probability theory — 12 , 8| . These ideas assume an earlier analysis 
of the evolution equation which, in the case of the Replicator-Diffusion model, has been 
made rigorous in [61]. A skillful, but heuristic, study can be found in 23[].This approach has 
two clear advantages: it allows for earlier heuristically derivations to be recast in a rigorous 
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^Neutral Kimura 

Trivial 
Generalized Kimura 

1/2 

Undefined \ Replicator 



1/2 



Figure 1: Fitnesses scales as \J>W(x, N, At) = 1 + (At)" ip^(x) and the relationship between population 
size and time-steps is given by TV -1 oc (A/)^. The graph above shows all possible macroscopic limits 
At — > for the solution of the Moran process. 



framework; it does not depend, necessarily, on the stochasticity of the underlying discrete 
model. A clear restriction is that it is a mean-field type model and only large populations 
can be treated. 

1.3. Outline 

In section |2j we briefly describe the discrete Moran process, with special emphasis on the 
invariants of the process. In section[3], we obtain the formal limit of the Moran process in the 
large population limit. The precise limit will depend on the different scalings assumption. 
In particular, for a certain choice, we obtain a maximal balance model that takes into 
account both genetic-drift and natural selection: the replicator-diffusion equation ffT5l) . In 
section HI we review some of the results obtained in [6( regarding the replicator-diffusion 
equation. As a particular application of the results, we point out how a generalization 
of the Kimura equation to the frequency-dependent case can be obtained through duality 
from ffTo^) . We also prove the convergence of the discrete model to the continuous one, 
in the appropriate limit. We also present some additional results on the connection of 
this model with game-theoretic issues on optimal fixation strategies. In section O we 
study another class of limit models, and show that they are equivalent to the Replicator 
equation. Additionally, we also study the fixation probabilities in the small diffusion limit. 



2. Discrete Dynamics 

We consider a fixed size population with two types of individuals, A and B say. We want 
to study the following evolutionary process: at a given time, we chose a newborn as copy 
of one of the individuals taken from the current population with probability proportional 
to its fitness. This newborn replaces a randomly chosen individual which is, therefore, 
killed. The process is then repeated after a timestep of At. This is the so-called Moran 



process 24] : see Figure [2] for an illustration. 
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(b) 



Figure 2: The Moran process: from a two-types population (a) we chose one to reproduce with probability 
proportional to the fitness, ans a second one (possibly the same) to be eliminated with equal probability (b), 
generating a new generation (c). 



For each type, we associate a fitness function depending on the type frequency: 

^(A)^(B) . [ Q) l] ^R+. 

Let P(x, t, N, At) be the probability that there are n = xN, with x e {0, ^, • • • , 1}, 
individuals of type A at time t, in a population of fixed size N. Let At be the time-step 
between each re-sampling of the population. At each step, the size of the type A can either 
increase by one, stay constant or decrease by one. If we denote the transitions probabilities 
by c + (x, N, At), c (x, N, At) and c_(x, N, At), respectively, we find that: 

x(l-x)^ A \x,N,At) 

C+lX ' ' ' xVW{x,N,At) + {l-x)y^{x,N,At) ' { ) 

x 2 ^(x,N, At) 

c {x,M,At) - x y (A)(XiNAt) - (1 _ x) ^(b )(X) ^ Af) ( 4 ) 



(1 -x) 2 ^ M \x,N, At) 

+ x y(A)(x,N,At) + (1 _ X )W{X, N, At) ' 

c (x N At) = x(l-x)^\x,N,At) 

From that, we may easily write an equation for the evolution of P: 

P(x,t + At,N,At) = c + (x-A^- 1 ,A^,At)P(x-A^- 1 ,t,iV,At) 

+c (x, N, At)P(x, t, N, At) 
+c_(x + N~\ N, At)P(x + N~\ t, N, At), 

x G {0, 1/N, . . . , 1} . (6) 

Let us define the relative fitness 

, . i ^ A \x,N, At) 

Also, let 

f(x) = x(l — x) and g(x; p) = 1 + (p — l)x . 
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Then it is a straightforward computation to verify that 

f(x)p(x,N,At) 



c+(x,N, At) 
Co{x,N, At) 
c-(x,N, At) 



g{x,p{x,N, At)) 

m 



i - 



g(x,p(x,N, At)) 
/(*) 



(l + p(x,N,At)) 



Let 



g{x,p(x,N,At)) ' 



1 



P(t)=[P(p,t,N,At),P[jj,t,N,At 



P(l,t,N,At) 



t 



Then we can rewrite fl6]) in matrix form as 

P(t + At) = MP(t), 
where M is a (N + 1) x (N + 1), tridiagonal matrix, with entries given by 



Ma = c ( — ,N,At J , i = 0,...,N 



M 



c, — ,iV,At z = JV-1, 



M 



— ,N,At 



i = 0,...,N- 1. 



Since iV, At) > 0, it is easy to see that M is a non-negative matrix. Since c_(x, N, At) + 
Co(x, A/", At) + c + (x, N, At) = 1, M is column stochastic. 

The long term evolution is described by the following theorem: 



Proposition 1. Let M be as above. Then 
1. 



(\ 1 - F 1 




lim M 







( 

\0 ^ ... Fjv_! 1/ 

where the F n satisfy 

F n = c+ iV, At) F n+1 + c (jj, N, At) F n + c_ JV, At) F n _ x , 
F = and Fjv = 1. 



(8) 



In particular, any stationary state will be concentrated at the endpoints. 
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2. If 1 denotes the vector (1,1,..., l)t ; F = (Fq, F±, . . . , Fn)' and if (•, •, ) denotes the 
usual inner product, then we have that 

(P(t),l) = <P(0),l) and (P(t),F) = (P(0),F). 

In particular, the I 1 -norm of a non-negative initial condition is preserved. 

Proof. For part 1, it is clear that from the definition of the transition matrix M that 
e = (1, 0, 0, ... , 0) and e N = (0, 0, . -j, 0, 1) are eigenvectors associated to the eigenvalue 
1. From the Gersgorin theorem (see [16( for instance), it is clear that all eigenvalues A are 
such that |A| < 1. For tridiagonal matrices with positive entries, the spectrum is real, i.e., 
o-(M) C R. Let us consider the matrix M := M N . Then, M 00 = M NN = 1, M (m)0 = 
~M iN = for i = 0, . . . , iV - 1 and M id > for i = 0, . . . , N and j = 1, . . . , N- 1. This last 
inequality represents the fact that it is always possible to go from any non-trivial initial 
condition to any other possible state in N time-steps. In particular, < Y2f=i ^-ji < 1 
fori= 1,...,N- 1. ^ 

We now prove that any stationary state (i.e. an eigenvector associated to an eigenvalue 
A = 1) of the matrix M must be a linear combination of eo and ejv- Therefore, this will 
be also true for M. Let P be a state with P = Pn = 0. Then 

AT-l JV-l N N-l N-l 

E ( M¥ )j = E = E ^ < E ■ 

Finally, MP ^ P. As for part 2, we first observe that, if a vector V satisfies IVr V = V, 
then we have that 

(P(t + At), V) = (MP(t),V) = (P(t),MtV) = (P(t),V). 

Hence (P(t), V) = (P(0), V). From the fact that M is column stochastic, we easily conclude 
that M^l = 1, and the first invariant follows. For the second invariant, we observe that 
Equation (jSJ) can be written in matrix notation as M'F = F, which concludes the proof. □ 

Remark 1. The invariants described in part 2, while natural at the discrete level, will play 
an important role in the correct specification of the continuous problem. 

The F n turns out to be the fixation probability of type A individuals, when the process 
starts in t = with n type A individuals. 

From the definitions of c*(x, N, At), * = — ,0,+, and p n = p (jj, N, At), we see that 
F n satisfies: 

p n F n +i — {p n + 1) F n + F n _i = , , s 

F = and Fn = 1. 1 ' 

Equation (Q can be solved exactly, with solution given by 

2^fc=i 1L=i Pi 

The expression given by ffTU]) does not appear to yield a simple formula in the general case. 
For constant fitness, see, however, [2j 
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3. Scalings and Thermodynamical Limits 

3.1. Preliminaries 

The aim of this section is to derive a continuous approximation, i.e., a PDE model that 
approximates the discrete process described in Section [2], in the large population limit. 

We define the probability density that at time t we have a fraction x e [0, 1] of type A 
individuals 

. P(x,t,N,At) , n 

V{x,t,N,At):= K ^, with x = -, n = 0,l,2,...,iV. 

Furthermore, we assume that in the thermodynamical limit, V(x, t, N, At) converges, in a 
sense to be made precise in section 4.2, to a function p(x,t) which is sufficiently smooth 
so that 

p(x±jj,t)= p(x, t) ± ±d x p(x, t) + ^d 2 x p(x, t) + 0(N~ 3 ) , (11) 
p{x,t + At) = p{x, t) + (At) d t p{x, t) + ((At) 2 ) . (12) 

3.2. Selection scalings 

From the weak selection hypothesis, equation (pQ), we write the relative fitness, defined 
by equation ([7|), asymptotically as 

p(x, N, At) « 1 + (At)" (^^(x) - ^ (B) (x)) + o ((At)") . (13) 

We also find up to order y (At)", defining y = iV -1 , 

c+(x-y,N,At) wi(l-i) [l + (l-x) (At)" (^ A )(x) - ^ (B) (x))] + y(l-2x) 
- y 2 - (At)"y [x(l - x) 2 d x (^ {A) (x) - ^ B \x)) -x(l-x) (^ {A) (x) - ^ B) (x)) 
+ (1 - 2x) (1 - ar) (^ (A) (a;) - ^ (B) (x))] , 



c (x, JV, At) « 1 - 2x(l - x) 



l-(At)"l^(^ A) (x)-^ (B) (x)) 



c_(x + y, iV, At) » z(l -x)[l-x (At)" (^ (A) (^) - ^ (B) 0))] + y(l - 2x) - y 2 
— (At)"y[x 2 (l — x)d x (ip( A \x) — ip^(x)) + x(l — x) (^)( x )-^( x )) 
+ [l - 2x)x (^ K \x) - ^ K \x))] , 

Joining these last three equation to the evolution equation (jBJ) and using assump- 
tions (dU) and flTJ), we find 

$p = -T^p;^ - *) (^ (A) (^) - ^ {B) (^)) P) + T^dl (*(1 - x)p) . 
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Now we assume N~ l = y = e (At) 11 , with e > 0, and take the limit At — > and, re-scaling 
time, we find all non-trivial possibilities according to the Figured], where for /i e (|>l]> 
v — 1 — n we find the diffusion equation 

dtp = \dl (x(l - x)p) ; (14) 
for jj, — v — |, we find the replicator- diffusion equation 

d t p = e -d 2 x (x(l - x)p) - <9 X (x{l - x) (^ A \x) - ^ (B) (x)) p) ; (15) 
and finally for v > fx = | we find the (partial differential version of the) replicator equation: 

d t p=-d x (x(l-x)(ij w (x)-tp m (x))p) . (16) 
The invariants in Proposition [T] become the following conservation laws: 

d f 1 

— / p(x, t) dx = 0, 
dt7o 



(17) 



d f 1 

— J 4>(x)p(x, t) dx = 0, (1^ 



o 



where 4> satisfies 

e 



-0"+(^(x)-^ B )(x))0' = O, 0(0) =0, 0(1) = 1. 

The conservation law (fT7|) applies to ffl4|) . ffT5|) and (T3~6]) • The conservation law ffl8|) applies 
to dH]) (with <p(x) = x), and to ffT5]) . 



4. The Replicator-Diffusion Equation 

The equation ffl5|) . without the conservations laws (fTT|) and ffl8|) . is a classical parabolic 



problem that is well posed, without imposing any boundary conditions [13|, |10j. However, 
it also true that any such a solution will vanish in the long time limit and, thus, f[T7|) and 
ffTS]) cannot be satisfied. 

Is possible to circumvent this difficulty, by using an appropriate notion of weak so- 
lution, that is somewhat different from the standard one. This has been done in detail 



in [6j, to where we refer the interested reader. In subsection 14.11 we briefly overview the 
results presented in jf|; in particular, we point out that a generalized form of the Kimura 
equation can be obtained as a dual of (Tl5l) and that, with this notion of solution, the 
forward/backward Kolmogorov pair is related through duality. In subsection 14.21 we show 
that the replicator-diffusion equation ffToT) is the right approximation of the discrete Moran 
process in the scaling u, — v — |, used through this entire section. Finally, in subsec- 
tion I4.3[ we deviate and, on assuming the usual relationship between fitness and payoff 
in game-theory, we study the connection of the Replicator-Diffusion fixation probabilities 
with strategies, and show that the strategies that maximize the fixation probability are 
exactly the same strategies that maximizes fitness, i.e., the Nash-equilibrium strategies. 
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4-1. The solution to the continuous problem 

The two conservation laws obtained in Proposition [1] generalizes for the continuous case 
at the form <9 t J p<pdx = 0, where <fi is any of the two linearly independent solutions of 

+(^ k \x)-^\x)) 0' = O . 

One solution is given by <pi(x) = 1 and the associated conservation law is simply the 
conservation of probability J pdx = 1. The second one is given by 



/; exp -f (^W ( x ") - ^ (B) (x")) dx" 



dx' 



Ms) = -i I -, ... T • (19) 



/o exp -f (^W {x") - {x")) dx" 



dx' 



Equation ffl5l) with the two conservation laws is well posed problem, if the appropriate 
notion of solution is used. See |6( for further details. Here, we only review the main result 
concerning the solution of this problem. 

Let BM + ([0 } 1]) denote the posit ive Radon measures in [0,1]. Then we have: 

Theorem 1. For a given p 1 e BA4 + ([0, 1]), there exists a unique (weak) solution p to 
Equation < f73|) . with p 6 L°° ([0, oo); BA4 + ([0, 1])) and such that p satisfies the conserva- 
tions laws d t Jq p(x, t)4>i(x)dx, i = 1,2. The solution can be written as 

p(x, t) = r(x, t) + a(t)8o + b(t)5i, 

where r G C°° C°°([0, 1])) is a classical (regular) solution to fTS\) without boundary 
conditions, and 8 y denotes the singular measure supported at y. We also have that a(t) 
and b(t), belong to C([0, oo)) fl C°°(1R + ). In particular, we have that 

p G C°°{R + , BM + {[0, 1])) n C°°{R + , C°°((0, 1)) . 

For large time, we have that lim^ 00 r(x,t) = ; uniformly, and that a(t) and b(t), the 
transient extinction and fixation probabilities, respectively, are monotonically increasing 
functions. Moreover, we have that 

lim p(t, •) = 7To[pVo + ^i[p%, 

t—>oo 

with respect to the Radon metric. Finally, the convergence rate is exponential. 

The precise values of tt and 7Ti (respectively, the final extinction and the fixation 
probabilities) are given by 7r [p 1 ] = 1 — ^[p 1 ] and iri[p l ] = p l (x)4>2(x)dx. 

Moreover, under this notion of solution, the adjoint equation to (|T5|) is given by 

d t f = x(l - x)d 2 J(x, t) + x(l - x) (^ A \x) - ^ M \x)) d x f(x, t), 
are (0,1), t>0, 

f(0,x) = f°(x). (20) 



Thus, through duality, equation ( 120]) generalizes the classical Kimura equation to the fre- 
quency selection case. In a particular case, this has been already presented in [sj. More 
recently, using a different approach, Champagnat et al. have also obtained a generalization 
of the Kimura equation to frequency dependent selection 
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4-2. Convergence of the discrete model to equation ( TI3)) in the limit of large populations 

In this subsection, we rigorously show that: take a general initial condition and evolve 
it through the replicator-diffusion equation ( fT5l) until time t and obtain the solution p CO nt; 
consider a discretization in intervals of size 1/N of the initial condition and evolve it 
through the Moran process -h = N 2 t times and obtain p N . Then the large population 
limit N — > oo of pn is, for any time t, pointwise, the solution of the continuous model p CO nt- 
This show that the continuous equation ffTol) is the right limit of the discrete process 
(Moran). In a sense, our approach can be described as inverse numerical analysis, as we 
look for a continuous model that approximates the discrete process, rather than the other 
way round, which is usually the point of view of numerical analysis. Then, it comes as no 
surprise that, we are able to prove convergence by using a classical theorem in numerical 
analysis. 

Theorem 2. Let p(x,t, N, At) be the solution of the finite population dynamics (of popula- 
tion N, time step At = 1/N 2 , i.e., /i = \), with initial conditions given by p°(x, N, At) = 
p°(x), x = 0, 1/N, 2/N, ■ ■ ■ ,1, for p° as in theorem^ Assume also the weak-selection 
limit, i.e, equation f[T5|) . with v — |. Let p cont (x,t) be the solution of the continuous model 
given by equation l[T5\) . with initial condition given by p°(x). If we write p™ for the i-th 
component of p(x,t, N, At) in the n-th iteration, we have, for any t* > 0, that 



lim p% 2 = p cont (x, t), x 6 [0, 1], t e [0, r 

N— >oo 

Proof. First, we consider the matrix M obtained from M by deleting the first and last rows 
and columns. Then, we observe that^the derivation of the thermodynamical limit shows 
that the discrete iteration given by M is consistent — in the approximation sense [33J - 
with Equation ffTBT) . without any boundary conditions, provided that we use (fT5|) . From the 
results of Proposition [TJ we know that the discrete iteration is stable, since cr(M) C (—1, 1). 
From Theorem [Tj, we see that the continuous problem without boundary conditions is well 
posed. In this case, we can then invoke the Lax-Ricthmyer equivalence theorem [33(] to 
guarantee that the discrete model converges to the continuum one, in the limit At, Ax — > 0, 
with At = (Ax) 2 . More precisely, the iteration defined by M converges to r(x,t), the 
smooth part of p CO nt(x,t); cf. 0]. 

Now returning to the iteration defined by M. In order to finish the proof, we only need 
to show that the first and last rows of M converge weakly to the appropriate Dirac masses. 
We shall do the computation for x = 0, the case x = 1 being similar. 

Since e G W N+1 converges weakly to 5q as N — * 00, we need only to show that it has 
the correct mass at each time t. 

For x = the iteration defined by M reads 

P(0, * + At, N, At) = P(0, t, N, At) + (^,t, N, At 

Thus, letting t = and solving the recursion, we have that 

1 m_1 / 1 \ 
P(0, mAt, N, At) = P(0, 0, N, At) + — ^ P ( -J At, N, At J 
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Dividing by 1/N, and taking the limit, we find, in a weak sense, that 



N 



lim Pq 



r(s, 0)ds + a 5o 



□ 



4-3. Strategy dominance 

For small populations it is know that the strategy that maximizes the fitness do not 



necessarily is the same as the strategy that maximizes the fixation probability [29] . In 



this section we identify types and strategies, as defined by game-theory [15|, |34|, and show 
that, for large populations, any strategy that maximizes the fitness will also maximize the 
fixation probability. For a study for finite populations similar to the one developed here, 
but that do not extend to infinite populations, see 26 . 

We say that a strategy A dominates a strategy B, A y B, if its fixation probability 
TTifp 1 ] > vr^fp 1 ] for any initial condition p 1 , where irf [p 1 ] = p 1 (x)xdx is the neutral 
fixation probability. This last value is obtained imposing = ip^ in equation (|19p and 
then 02 (%) = x. It is clear that the relation y is reflexive, anti- symmetric, and transitive; 
i.e., y is a partial order in the set of all possible strategies. This cannot be made a total 
order as there are games, like the Hawk-and-Dove game, where no pure strategy dominates 
the other. 

First note that A y B if and only if tti^o] > irf [5 x o] for x Q G [0, 1]. Now, let us suppose 
that A is the Nash equilibrium strategy. Then, for any strategy B and for any particular 
composition of the population, ip( A \x) — ^ M \x) > 0. This implies that the function 
J* (ip( A \x') — ip( B \x')^ dx' is a non-decreasing function of x E [0, 1] and consequently 



x 



exp 



dx' 



> 



r 1 


2 r' 


/ ex p 




Jo 


£ Jo 



dx' . 



(^ A \x")-^ M \x"))dx" 
(^ A \x")-^ M) (x"))dx' 
In particular, ^{x ) > x , Vx° G [0, 1] and finally 



7Ti [5 x o] — 7r} N [8 x o] = I 5 x o(x) (4>2(x) — x) dx = 4>2(x°) — x° > . 
'o 



5. The Replicator Dynamics 

In subsection 15.11 we show how the Replicator Dynamics can be obtained from the 
discrete Moran process, when the scalings are chosen such that the (fT6j) is the correct 
limit. In subsection 15.21 we study the fixation probabilities in the small diffusion limit. 
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5.1. The replicator equation as a large population limit 

The replicator dynamics models the continuous time evolution of the fraction of a given 
type of individuals in a infinite well-mixed population framework. In its simplest form, it 
reads as follows: 

X = X(1-X)(?/; (A) (X)-^ (B) (X)) . (21) 

We shall now see that, when the scalings are such that ffl6|) is the correct limit, we actually 
obtain the Replicator dynamics. We also want to point out that the Replicator dynamics 
can also be rigorously obtained as the small diffusion limit, or, equivalently, the limit 
of strong selection and intermediate times, of the replicator-diffusion equation ( fl5l) . For 
further details, see [51. 



Let 

u{x,t) = x(l ~ x)(ip {A) (x) - ip m (x))p(x,t) 
Then, the equation ( fl6l) can be written as 

d t u = x(l - x)(</> {A) (x) - ^ M \x))d x u 



Its characteristics [18( are given by 

^ = x{l-x) U^{x)-^ v \x)) 
as 

The projected characteristics in the x x t plane are, therefore, given by 

which is exactly the equation f[2"Tj) . 

The solution, for smooth initial condition p°(x), can be written as 

, A o ((h , , ^- t (x)(l-$- f (x)))(^ A )($_ f (x))-^)($_ f (x))) 

p( X , t) = P x^xwm^wm ■ (22) 

where is the flow map of (I2ip . and the fraction is interpreted as being one at every 

equilibrium point of (|2~TT) . 

Thus, equation (FTO]) can be seen as a Eulerian representation of u, while ( 12~T|) would be 
a Lagrangian representation. From u to p there is a bijective map. Thus, ffT6]) and fT2~T]) 
are equivalent. 

The (121]) can be SGG clS db Lagrangian representation of u, once the initial probability 
distribution is given, while ffTBl) can be seen as an Eulerian representation of the same 
quantity. 

Remark 2. A similar version of Theorem^ is valid for the convergence of the Moran 
process to (1161) . in the case of smooth initial conditions. The proof is very similar and, 
hence, is omitted. 
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Fixation probability 



A 










E=0 










/'e>0 










1 

— s- 



x* Initial presence 



Figure 3: Fixation probability 4>2(x) versus the initial presence x (i.e., initial condition is given by p l — 
S x The solid line shows the (discontinuous) lim E ^o </>2 , when there is a unique interior point such that 
i/jW(z*) = ^)( x *). This is the case of games like the stag-hunt or the hawk-and-dove. The dashed lines 
indicate 4>2 for positive e. Compare with figure 1 in (35|. 



5.2. Fixation probability and the small diffusion limit of equation J73| ) 

In this subsection we obtain the small diffusion limit e — > of the fixation probability fa 
given by equation (TT9]) . Let us consider an initial condition given by p l = 5 x o and then the 
fixation probability is given by fa (x°). Now, we obtain the limit e — > from equation (fl9l) . 
The only important terms when e — > come from the minima of (xp^ {%') — ip B ^ (x')j dx', 
see 



14i | . If the minimum x* is unique then 

\imfa(x) = X[o,x]( x *) 



where xi is the characteristic function of the interval I. In particular, if ^ > > tp™> then 
x* = and lim £ _> (f) 2 {x) = 1, for x G (0, 1), meaning that the type A will be fixed with 
probability 1, if it is initially present, even in a very small quantity. If there is a unique 
interior point, x*, such that both pay-off are the same, type A will be fixed if its initial 
presence is superior to x* and extinct otherwise. Note in figure [3] that the ^-diffusion is a 
regularizing effect in the fixation probability. Compare with Figure 1 in 35 1. See also j^]. 



If there are more than one minimum for (ip(^ [x') — tp M \x')) dx', the fixation proba- 
bility in the limit e — > is the average of the function X[o,x] over all minima. For example 
if /; (^ (A) (x') - ^ B) (x')) dx' = x(l - x), then linw, fa(x) = \ (x[o,*](0) + X[o,x](l)) = \ 
for x G (0, 1). In the neutral case, lim e ^ fa(x) = J* X[o,x](x')dx' = J* dx' = x. 



2 ' 



6. Discussion 

The results presented here can be seen as an alternative approach both to mean-field 
approximations usually made through diffusion theory [12J and to more recent probabilistic 
approaches [8|. 



14 



While the analysis presented was focused on the Moran process, the paradigm pre- 
sented here can accommodate more general discrete models. The formal derivation of the 
continuous model — in a self-consistent manner — as a first step, should be seen as a feature, 
since it allows for the possibility that a heuristically derived approximation can be made 
rigorous. It is also less dependent on the stochasticity of the subjacent discrete model, 
which can be an advantage in many situations. 



It also shows that diverse models as evolutionary game dynamics [15l . |30T | and neutral 



genetic variability [19|, [12| can be tighten from the discrete process. As quoted in [36 
while Darwinian motto is the survival of the fittest, Kimura work on neutral selection in 
genetics was described by himself as survival of the luckiest. Here we show that both can 
be accommodated — see also 0] . Moreover, we show that from a unique discrete process, 
we can have both effects, and recover either Darwinian, pure Stochastic, or a balanced 
evolutionary dynamics. A more detailed study of the interconnection between Darwinian 
and Stochastic evolution as a multiscale process is being carried by [5j. This also suggests 
that an important issue in Mathematical Evolutionary Dynamics is to understand the 
interplay between the effects and time scales in the diverse biological phenomena linked to 
the evolutionary process. 
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